ON THE TWO-PHASE NAVIER-STOKES EQUATIONS 
WITH SURFACE TENSION 
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Abstract. The two-phase free boundary problem for the Navier-Stokes 
system is considered in a situation where the initial interface is close to 
a halfplane. By means of Lp-maximal regularity of the underlying linear 
problem we show local well-posedness of the problem, and prove that 
the solution, in particular the interface, becomes instantaneously real 
analytic. 
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1. Introduction and Main Results 

In this paper we consider a free boundary problem that describes the motion of 
two viscous incompressible capillary Newtonian fluids. The fluids are separated by 
an interface that is unknown and has to be determined as part of the problem. 
Let fii(O) C ]R"+^ (n > 1) be a region occupied by a viscous incompressible fluid, 
fluidi, and let fl2{0) be the complement of the closure of ili(O) in corre- 
sponding to the region occupied by a second incompressible viscous fluid, fluid2. 
We assume that the two fluids are immiscible. Let Fq be the hypersurface that 
bounds ili(O) (and hence also il2(0)) and let r{t) denote the position of Fq at 
time t. Thus, T{t) is a sharp interface which separates the fluids occupying the 
regions fli{t) and il,2{t), respectively, where (t) := R"+^ \ ili(t). We denote the 
normal field on F(t), pointing from ni{t) into il.2{t), by v{t,-). Moreover, we de- 
note by V{t, ■) and ^(t, •) the normal velocity and the mean curvature of T{t) with 
respect to v{t, •), respectively. Here the curvature k{x, t) is assumed to be negative 
when rii(i) is convex in a neighborhood of a; e F(i). The motion of the fluids is 
governed by the following system of equations for « — 1, 2 ; 

' -I- (ii|V)ii) - ^iAu + Vg = in Oi(t) 

div u = in 17, (t) 

— \S{u,q)v\=aKV on T{t) 

M=0 on V{t) (1.1) 

V = {u\v) on T{t) 

m(0) = Uq in VL,{0) 

r(o) = Fo. 
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Here, S = S{u, q) is the stress tensor defined by 

S{u,q) ^ fi^{Vu+ {Vu^) ~ ql in n,{t), 

and 

denotes the jump of the quantity v, defined on the respective domains Qi{t), across 
the interface T{t). 

Given are the initial position Tq of the interface, and the initial velocity 

The unknowns are the velocity field u{t, •) : ^l{t) K"+^, the pressure field q{t, •) : 
n{t) R, and the free boundary T{t), where fl{t) := fJi(t) U fl2it). 
The constants pi > and pi > denote the densities and the viscosities of the 
respective fluids, and the constant cr stands for the surface tension. Hence the 
material parameters pi and pi depend on the phase i, but otherwise are assumed 
to be constant. System comprises the two-phase Navier- Stokes equations with 
surface tension. The first equation in (|l.ip refiects balance of momentum, while the 
second expresses the fact that both fluids are incompressible. If surface tension is 
neglected, the boundary condition on T{t) would be the equality of stress on the 
two sides of the surface. The effect of surface tension introduces a discontinuity 
in the normal component of proportional to the mean curvature of T{t). 

The forth equation stipulates that the velocities are continuous across T{t). Finally, 
the fifth equation, called the kinematic boundary condition, expresses that fluid 
particles cannot cross T{t). 
In order to economize our notation, we set 

P = PiXniit) + P2Xn2(t), P = PiXniit) + P2Xn2{t), 

where xd denotes the indicator function of a set D. With this convention, system 
(jl.ip can be recast as 

' p{dtu + {u\V)u) 



(1.2) 



— pAu + Vg 
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in 


n{t) 


div u 
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In this publication we consider the case where To is a graph over M" given by a 
function Hq- We then set rii(O) — {{x, y) £ M" x M. : y < ho{x)} and consequently, 
^^2(0) = {{x,y) G M" X R : y > ho{x)}. Our main result on existence, uniqueness, 
and regularity of solutions then reads as follows. 

Theorem 1.1. Suppose p > n + 3. Then given to > 0, there exists So = £o{to) > 
such that for any initial values 

{uo,ho) e M^p2-2/p(^^^^n+l) X VFp3-2/P(R"), 
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satisfying the compatibility conditions 

l^j.D{uo)vo - K^o\D{uo)voyaj = 0, div mq = on fta, (uoj = 0, 
with D(uq) :— (Vuq + (Vuo)'''), and the smallness condition 

problem (|1.2p admits a classical solution F) on (0, to)- T'^e solution is unique 
in the function class described in Theorem \ 6.3[ In addition, T{t) is a graph over 
given by a function h{t), A4 — Ute(o to) ({^1 ^ ^(t)) is a real analytic manifold, 
and with 

O ^{{t,x,y): te (0, to), a; e M", y ^ h{t, x)}, 
the function (u, q) : O ^ R"+^ is real analytic. 

Remarks 1.2. (a) Theorem 1.1 shows that solutions immediately regularize and 
become analytic in space and time. If one thinks of the situation of oil in contact 
with water, this result seems plausible, as capillary forces tend to smooth out corners 
in the interface separating the two different fluids. 

(b) More precise statements for a transformed version of problem (II. 2p will be given 
in Section 6. Due to the restriction p > rt + 3, we shall show that 

h^C{J■BUC'^{W'))C^C^{J■BUC^{W')) (1.3) 

where J — [0, to]- In particular, the normal of ^i(t), the normal velocity of T{t), 
and the mean curvature of V{t) are well-defined and continuous, so that ()1.2|) makes 
sense pointwise. For u and q we obtain 

u{t,-) e BUC\n{t),R"+^) fort e J, ue BUC{J X IR"+\R"+1), 

g(t,-) e [/C(f^(t)) fort e J\{0}. ^^'^^ 

In addition, the solution (u, q, h) depends continuously on the initial values (wq, ho). 
Also interesting is the fact that the surface pressure jump will turn out to be real 
analytic as well. 

(c) It is possible to relax the assumption p > n + 3. In fact, p > (n + 3)/2 can be 
shown to be sufficient. However, to keep the arguments as simple as possible, here 
we impose the stronger condition p > n -I- 3. 

(d) If gravity acts on the fluids then the condition on the free boundary is to be 
replaced by 

-lS{u,q)}v^aHv + ^lp\yv on F(t), (1.5) 

where y denotes the vertical component of a generic point on F(t), and where 7 > 
is the gravity acceleration. It is not very difficult to verify that our approach also 
covers this case, yielding a solution having the same regularity properties as stated 
in the theorem above, provided p2 < Pi, i-S- the heavier fluid lies beneth the lighter 
one. Indeed, an analysis of our proof shows that we only need to replace the symbol 
s(A,t) introduced in (|5.9p by 

s(A, r) = A + crTfc(z) - ^^k{z). 

T 

It satisfies the same estimates as in ()5.10p in case that /O2 < Pi - 

(e) We mention that our results also cover the one phase Navier-Stokes equations 
with surface tension (II. 6p . 

(f) The solutions we obtain exist on an interval (0,to) with tg > arbitrary, but 
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fixed, provided the initial data are sufficiently small. It can be shown that prob- 
lem (|1.2p also admits unique local solutions that enjoy the same regularity properties 
as above, provided supj.gjj |V/io| is sufficiently small in relation to the horizontal 
component of mo. In this case, no other smallness conditions on the data are re- 
quired. The proof of this result is considerably more involved, and the analysis re- 
quires delicate estimates for the nonlinear terms. Additionally, we need a modified 
version of Theorem 5.1 in order to dominate some of the nonlinear terms by linear 
ones. The proof of this modification will involve introducing a countable partition 
of unity and then establishing commutator estimates for certain pseudo-differential 
operators. Since this paper is already rather long, we refrain from including a proof 
of this result here. It will be contained in the forthcoming paper |32j . 

Let us now discuss and contrast our results with results previously obtained by other 
researchers. In case 5^2 (i) = one obtains the one-phase Navier-Stokes equations 
with surface tension 

( p{dtu + iu\V)u) 



(1.6) 



Equations (|1.6p describe the motion of an isolated liquid which moves due to cap- 
illary forces acting on the free boundary. 

Problem (|1.6p has received wide attention in the last two decades or so. Existence 
and uniqueness of solutions for a > 0, as well as for cr = 0, in case that n{0) is 
bounded (corresponding to an isolated fluid drop) has been extensively studied in 
a long series of papers by Solonnikov, see for instance [SH]"!!!] and [57] for the 
case CT > 0. Solonnikov proves existence and uniqueness results in various function 
spaces, including anisotropic Holder and Sobolev-Slobodetskii spaces. Moreover, it 
is shown in |40j that if f^o is sufficiently close to a ball and the initial velocity mo 
is sufficiently small, then the solution exists globally, and converges to a uniform 
rigid rotation of the liquid about a certain axis which is moving uniformly with 
a constant speed, see also [28 . More recently, local existence and uniqueness of 
solutions for (|1.6[) (in case that is a bounded domain, a perturbed infinite layer, 
or a perturbed half-space) in anisotropic Sobolev spaces W^^^ with 2 < p < oo and 
n < q < oo has been established by Shibata and Shimizu in [37l [38] . For results 
concerning p.6|) with cr = we refer to the recent contributions [35l [SG] and the 
references therein. 

The motion of a layer of viscous, incompressible fluid in an ocean of inflnite extent, 
bounded below by a solid surface and above by a free surface which includes the 
effects of surface tension and gravity (in which case JIq is a strip, bounded above 
by Fq and below by a fixed surface F;,) is considered by AUain [T], Beale [5], Beale 
and Nishida [7] , Tani [47] , and by Tani and Tanaka [48] . If the initial state and the 
initial velocity are close to equilibrium, global existence of solutions is proved in [6] 
for cr > 0, and in [48] for cr > 0, and the asymptotic decay rate for t — > oo is studied 
in [I]. 
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Results concerning the two-phase problem (|1.2p are more recent. Existence and 
uniqueness of local solutions is studied in [lOl [TTl [121 US- rnore detail, Den- 
siova [11] establishes existence and uniqueness of solutions (of the transformed 
problem in Lagrangian coordinates) with v G W2^^^ for r G (5/2,3) in case that 
one of the domains is bounded. Tanaka [46j considers the two-phase Navier-Stokes 
equations with thermo-capillary convection in bounded domains, and he obtains 
existence and uniqueness of solutions with {v,9) e W"^"^^"^ for r e (7/2,4), with d 
denoting the temperature. 

The approach used by Solonnikov, and also in [10]-[12|, [35l |Ml EH H Si HZl EH] , 

relies on a formulation of systems ()1.2p and ()1.6p in Lagrangian coordinates. In this 
formulation one obtains a transformed problem for the velocity and the pressure 
on a fixed domain, where the free boundary does not occur explicitly. The free 
boundary can then be obtained by 



where v is the velocity field in Lagrangian coordinates. It is not clear whether this 
formulation allows one to obtain smoothing results for the free boundary, as the 
regularity of r(t) seems to be restricted by the regularity of Fq. To the best of 
our knowledge, the regularity of the free boundary for the Navier-Stokes equations 
with surface tension ()1.1|) or (|1.6|) has not been addressed in the literature before, 
with the notable exception of [S]. Beale considers the ocean problem with f2(i) = 
{{x,y) € X M : —h{x) < y < h{t, x)} and he shows by a boot-strapping argument 
that solutions are C'^ for any given fixed fc e N, where the size of the initial data 
must be adjusted in dependence of k. As in our case, his approach does not rely on 
a formulation in Lagrangian coordinates. 

In order to prove our main result we transform problem (|1.2p into a problem on 
a fixed domain. The transformation is expressed in terms of the unknown height 
function h describing the free boundary. Our analysis proceeds with studying solv- 
ability properties of some associated linear problems. It is important to point out 
that we succeed in establishing optimal solvability results (also referred to as as 
maximal regularity), see Theorem 13.11 Proposition [3?3l Theorem 14.11 CoroUarv 14.21 
and Theorem 15. II In other words, we show that the linear problems define an iso- 
morphism between properly chosen function spaces. This property, in turn, allows 
us to resort to the implicit function theorem to establish the analyticity of solu- 
tions to the nonlinear problem, as will be pointed out below. All our results for the 
associated linear problems mentioned above seem to be new, as they give sufficient 
as well as necessary conditions for solvability. Our analysis is greatly facilitated by 
studying the Dirichlet-to-Neumann operator for the Stokes equations, see Section 4. 
It is interesting, and maybe even surprising, to observe the mapping properties of 
this operator, see Theorem 4.1. Our approach for establishing solvability results re- 
lies on the powerful theory of maximal regularity, in particular on the iJ°°-calculus 
for sectorial operators, the Dore-Venni theorem, and the Kalton-Weis theorem, see 
for instance [fl |M US 130] ■ 

Based on the linear estimates we can solve the nonlinear problem by the con- 
traction mapping principle. Analyticity of the solution is obtained in a rather short 
and elegant way by the implicit function theorem in conjunction with a scaling ar- 
gument, relying on an idea that goes back to Angenent [31 [4] and Masuda [26]; see 
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also [171 UHl EQ] . More precisely, by introducing parameters which represent scal- 
ing in time, and translation in space, the implicit function theorem yields analytic 
dependence of the solution of a parameter dependent-problem on the parameters, 
and this can be translated to a smoothness result in space and time for the original 
problem. 

The plan for this paper is as follows. Section 2 contains the transformation of 
the problem to a half-space and the determination of the proper underlying linear 
problem. In Sections 3, 4 and 5 we study this linearization and prove in particular 
the crucial maximal regularity results in an Lp-sctting. Section 6 is then devoted 
to the nonlinear problem and contains the proof of our main result. 



2. Reduction to a Flat Interface 

In this section we first transform the free boundary problem ()1.2p to a fixed 
domain, and we then introduce some function spaces that will be used throughout 
the paper. Suppose that r(^) is a graph over R", parametrized as 

T{t) = {(x, h{t, x)) : xe M"}, t e J, 

with n2{t) lying "above" r(t), i.e. ^2(0 = {{x,y) e W x R : y > h{t,x)} for 
t E J :— [0,a]. Reduction from deformed into true halfspaces is achieved by means 
of the transformations 

ui{t,x,h{t,x) +y) 
v{t,x,y)= ■■ 

_ Un{t,x,h(t,x) +y) 
w{t,x,y) = Un+iit,x,h{t,x) +y), 
T^{t, X, y) = q{t, X, h{t, x) + y), 

where t E J, x E M", y G R, y 7^ 0. Since for j,k ^ 1, . . . , n we have 

djUk = djVk - djhdyVk, dn+iUk = dyVk, 
djUn+i = djW - djhdyW, 9„+im„+i = dyW, 
djQ = dj-K - djhdyTT, dn+iq = dyir, 
dtUk = dtVk - dthdyVk, 9fU„+i = dtw - dfhdyW, 



(2.1) 



and 



Auk = A^vk - 2{Vh\V^)dyVk + (1 + \Vh\^)d'^Vk - AhdyVk, 
Aun+i = A^w ~ 2{Vh\V^)dyW + (1 + \Vh\'^)dyW - AhdyW, 



we obtain from (|1.2p the following quasilinear system with initial conditions 

pdtv - ^iA.^v - ^id^v + V.^7r ^ Fy{v,w,'K,h) in (0, cxj) x R"+^ 
pdtw — ^A,j,w — ^dyW + dyiT — Fyj{v,w,h) in (0, 00) x R"^"^ 
AiY^v + dyW Fd{v,h) in (0, 00) x R"+i 



(2.2) 



■"(O, X, y) = va{x, y), w(0, x, y) = ^0(2;, y) 



where ~ {{x, y) e M" x M ; y 7^ 0}. Here and in the sequel, Vh and Ah always 

denote the gradient and the Laplacian of h with respect to x £ R". Note that p 
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and iJL in general have jumps at y = 0, i.e. p — p2 foi y > 0, p = pi for ?/ < 0, and 
similarly for fi. The nonlinearities are given by 

Fy{v, w, n, h) = ^i{~2{Vh\V:,)dyV + |V/ip9^i; - l\hdyv} + dyirVh 

+ p{—{v\\/ x)v + (V/i|u)i9yU — wdyv} + pdfhdyV, 

F^{v, w, h) = ^{-2(V/i|V,;)ayW + \Vh\^dlw - l^hdyw} (2.3) 

+ p{ — {v\\/ x)w + (V/i|w)9yU' — W(9j;w} + pdthdyW, 

Fd{v,h) = {\/h\dyV). 

Note that these functions are polynomials in the derivatives of {v,w,T:^h), hence 
analytic, and linear with respect to second derivatives, with coefficients of first 
order. This exhibits the quasilinear character of the problem. 

To obtain the transformed interface conditions we observe that the outer normal 
V of ^i{t) is given by 



Vit, X) = 



1 



where, as above, Wh{t, x) denotes the gradient vector of h with respect to a; S 
The normal velocity V of r(-) is 



V{t,x) =dthit,x)/.yi + \Vh{t,x)\'^. 

The kinematic condition V = on r(-) now reads as 

dth-jw = H{v,h), H{v,h) -.^ ~{jv\Vh). (2.4) 

Here {'jw){x) :— w{x,0) denotes the trace of the function w : M"+^ R and, 
correspondingly, is the trace of v : — > M". Since u is continuous across 
r(<), 7w and "fw are unambiguously defined. It is also noteworthy to observe that 
the tangential derivatives of v and w are continuous across M". The curvature of 
T{t) is given by 

see for instance equation (24) in [8, Appendix], with 

_ \Vh\^Ah iVh\V^hVh) 

' (i + v/rTW)v/rTw (i + iv/^P)'/'' 

where V^/i denotes the Hessian matrix of all second order derivatives of h. The 
components of 2?(w, w, h), the transformed version of the deformation tensor D{u) = 
(Vu + (Vu)^), are given by 

Vij (w, w, h) — diVj + djVi — {dihdyVj + djhdyVi), 

'Dn+ijivTW, h) = I?j,„+i («,!(;, h) = dyVj + djW - djhdyW, (2.6) 

2?»i+i,ri+i(w, w, /i) = 2dyW, 

for i,j — 1, . . . ,n, where Sij denotes the Kronecker symbol. For the jumps of the 
components of the deformation tensor this yields 

lfJ,V^j{v, w, h)j = l^l{^^VJ + ^jv^)j - ^^hl^i^yVJj - djhlfidyVij, 

lp.'Dn+ij{v,w,h)j = lfi'D.i,n+i{v,w,h)j = lfJ.djwj + llJ,dyVjj - djhlfidywj, 

lflT>n+l,n+l{v,W,h)j = 2[fidyWj. 
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Therefore, the jump condition for the normal stress at the interface yields the 
following boundary conditions: 



(2.7) 



-llJ.dyVj - IfiW^w] = G^{v,w, |7r],/i), 
-2lfidywj + Itt] - aAh ^ G^(v, w, h), 

where the nonlinearities (Gu, G^) have the form 

G,{v,w, liTlh)^-l^i{^^v + {^^vY)\-^h+\^h\''l^idyvl + (V/i| l^^^yv\)Wh 

- l^idywlVh + {H - a{Ah - G^{h))}Vh, (2.8) 
G^{v,w,h)^-{Vh\ l^l^yv\) - {Vh\ ItiVM) + \'^hWtidywl - aG^ih). 

We note that G — (G^,G„) is anaytic in (w,w, |7r],ft,). Moreover, G is linear in 
{v,w, [tt]), and in the second derivatives of h. Thus the boundary conditions are 
quasilinear as well. 

Summarizing, we arrive at the following problem for u — {v, w), tt, and h: 



pdtu — /iAu + Vtt 
div u 

-2yidywj + Itt] - aAh 

dfh — 'yw 
u{0) = 1*0, h{0) 



F{u, IT, h) 
Fd[u,h) 
G,(w, M,h) 
Gw{u, h) 


H{u,h) 
ho, 



m 
in 
on 
on 
on 
on 



pri+1 



(2.9) 



for t > 0. This is problem (|1.2p transformed to the half-spaces MJ^ 
R" xR:±y> 0}. 



Before studying solvability results for problem (|2.9p let us first introduce suitable 
function spaces. Let O C R™ be open and X be an arbitrary Banach space. By 
Lp{n;X) and Hp{il;X), for 1 < p < oo, s S K, we denote the X-valued Lebegue 
and the Bessel potential spaces of order s, respectively. We will also frequently make 
use of the fractional Sobolev-Slobodeckij spaces W^(S1; X), I < p < oo, s e M \ Z, 
with norm 



hWw^in-x) = \\9\\wl^}^n-x) + J2 (11 

\a\ = [s] 



1/p 

dxdy], (2.10) 



|9"g(x)~9"ff(y)||^ 

where [s] denotes the largest integer smaller than s. Let a € (0,oo] and J = [0,a]. 
We set 

r {g e W;{J;X) : .g(0) = .g'(0) - . . . = g^^HO) = 0}, 

oW;iJ;X):^ I if fc+i<.s<fc + l + i, /ceNU{0}, 

[ %^(J;X), if s<i. 

The spaces oHp(J;X) are defined analogously. Here we remind that — Wp for 
fc G Z and 1 < p < CX3, and that = for s e R \ Z. 

For C M™ open and 1 < p < oo, the homogeneous Sobolev spaces Hp{n) of 
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order 1 are defined as 

Him {{g e Li,ioc(f^) : ||V5|U^(o) < || • 

\\9\\Hi,r-=it\\9Mlmf^'- ^'-''^ 

Then Hp{fl) is a Banach space, provided we factor out the constant functions and 
equip the resuhing space with the corresponding quotient norm, see for instance \21\ 
Lemma II. 5.1]. We will in the sequel always consider the quotient space topology 
without change of notation. In case that f2 is locally Lipschitz, it is known that 
Hp{Q) C i?piQ(.(r2), see 'Qli Remark II. 5.1], and consequently, any function in 
Hp{n) has a well-defined trace on dH,. 

For s G M and 1 < p < oo we also consider the homogeneous Bessel-potential 
spaces i/p(M") of order s, defined by 

H;{R") {{g e 5'(R") : Pg G Lp(K")}, || • ||^.(«„)), 

(2.12) 

li5lli?|(R") — Wi'gh^iR"), 

where iS'(M") denotes the space of all tempered distributions, and I" is the Riesz 
potential given by 

Pg := i-Ay/^g .^-1(1^1^.^5), g G 5'(R"). 

By factoring out all polynomials, iJp(M") becomes a Banach space with the natural 
quotient norm. For s G IR\Z, the homogeneous Sobolev-Slobodeckij spaces Wp{M.") 
of fractional order can be obtained by real interpolation as 

W;{R^) := (ijp^(R"), iip^+i(R")),_fc,p, k<s<k + l, 

where (•, ■)e,p is the real interpolation method. It follows that 

r G Isom(ff^+'(R"),i7^(R")) n Isom(M7p*+''(M"), ^^^^(R")), s,i G R, (2.13) 

with = for fc G Z. We refer to Section 6.3] and [50l Section 5] for 
more information on homogeneous functions spaces. In particular, it follows from 
parts (ii) and (iii) in [5D1 Theorem 5.2.3.1] that the definitions (|2.1ip and (|2.12p are 
consistent if = R", s = 1, and 1 < p < oo. We note in passing that 

define equivalent norms on W^p(R") for < s < 1, where P(-) denotes the Poisson 
semigroup, see [501 Theorem 5.2.3.2 and Remark 5.2.3.4]. Moreover, 

7± G £(Typi(M!^+i),iypi-i/P(R")), (2.15) 

where 7± denotes the trace operators, see for instance [211 Theorem II. 8. 2]. 
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3. The Linearized Two-Phase Stokes Problem 
In this section we consider the hnear two-phase (inhomogeneous) Stokes problem 



pdtu — /iAm + Vtt 


= / 


in 




div u 


= fd 


in 


R"+i 




= Qv 


on 


R" 




= 9w 


on 


R" 


M 


= 


on 


R" 


u(0) 


= "0 


in 


R"+i 



Here the initial value uq as well as the inhomogeneities {f ^ fd, gvi 9w) are given. We 
want to establish maximal regularity for this problem in the framework of Lp-spaces. 
Thus we are interested in solutions {u, tt) in the class 

u G HliJ- Lp(R"+\ M"+i)) n Lp( J; i/2(R"+i^ r"+1))^ ^ ^ l^^j. iyi(M"+i)). 

We remind here that J = [0, a] and R"+i = {{x,y) e R" x M : y ^ 0}. If (u,7r) 
is a solution of (I3.1|l in this class we necessarily have / € Lp( J; Lp(R"+^)), and 
additionally uq G VKp"^/''(M"+\ R"+1) by trace theory. Moreover, 

as the operator div maps ip(]R"+^) onto _ffp"^(M"+^). Taking traces at the interface 
y ^ Q results in g„ e Wp^/^"^/^^( J; Lp(R", R")) n Lp( J; iyp~^/^^(R", R")), and 
g„ e Lp(J;iyp~^/^(R")). If, in addition, 

then (7u) shares this regularity. 

The main result of this section states the converse of these assertions, i.e. maximal 
Lp-regularity for p.ip . 

Theorem 3.1. Let 1 < p < oo he fixed, p ^ 3/2,3, and assume that pj and p,j are 
positive constants for j — 1,2, and set J = [0, a]. Then the Stokes problem p.ip 
admits a unique solution (u, tt) with regularity 

i?p\j;ip(M"+\R"+i))nip(J;i?p2(M"+\M"+^)), tt e Lp{J; H^{M.''+^)), 

if and only if the data (/, fdi dvi g-w, uq) satisfy the following regularity and compat- 
ibility conditions: 

(a) /eLp(J;Lp(R"+i,M"+i)), 

(b) fd e i7pi(J;ijp-H»"+')) n Lp(J;i/pi(R"+i)), 

(c) g, G iyp/'"'/'P(J;ip(R",R"))nLp(J;M^p'~'/P(R",M")), 
5„ eLp(J;Wp'-^/''(R")), 

(d) Mo G Typ~^/''(M"+i,R"+i), 

(e) div-uo = fd(0) m R"+i and |uo] = on R" i/p > 3/2, 

(f) -l^I^yVoj - I/iV,Wo] = 9viO) on M" if p > 3. 
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In addition, [tt] G Wp^^ ^/^^( J; ip(K")) n Lp( J; VKp ^/^(R")) if and only if 

The solution map [(/, /d, 5t), 5u), 5h, uq, /iq) '-^ (ujTt)] is continuous between the cor- 
responding spaces. 

Proof. The basic idea of the proof is to reduce system ()3.1|) to the case where 
(/, /d,uo) = (0,0,0) and .gi;(0) = 0, and then to solve the resuhing problem by 
means of the Dirichlet-to-Neumann operator for the Stokes problem. We can achieve 
this goal in four steps, as follows. 

Step 1. For given data (/, 5«,mo) subject to the conditions of the theorem we first 
solve the parabolic problem without pressure and divergence, i.e. we solve 



pdtu — fiAu 




in 






= gv 


on 




-ni^dywj 


= 9w 


on 


M" 




= 


on 


M" 


u{0) 


= Uq 


in 





Here we set g^, = — 2e~^"*|//9yWo] with Z3„ := —A in Lp(M"). The function 
has the same regularity as g^, and the necessary compatibility conditions are sat- 
isfied. By reflection of the {y < 0}-part of this problem to the upper halfplane, we 
obtain a parabolic system on a halfspace with boundary conditions satisfying the 
Lopatinskii-Shapiro conditions. Therefore, the theory of parabolic boundary value 
problems yields a unique solution ui for (|3.2p with regularity 

ui e ifi(J;ip(M"+\R"+i))nLp(J;i/2(M"+i,M«+i)). 

We refer to Denk, Hieber and Priiss [13l[14] for this. 
Step 2. In this step we solve the Stokes equations 

pdtU - pAu + Vtt = in M"+^ 

div M = /d - div Ml in M"+^ (3.3) 

m(0) = in M"+i 

where ui is the solution obtained in Step 1. It follows from assumption (e) that 
system (13.31) satisfies the compatibility condition divu(O) = /d(0) — divui(O) — 0. 
We remind that p = y02XR^+i +PiXm"+i and p, = ^2Xr^+i +MiXr"+i- Concentrating 
on the upper halfplane, we extend the function (/^ — divui) evenly in y to all of 
R"+^ and solve the Stokes problem with coefficients p2, P2 in the whole space, see [9l 
Theorem 5.1]. This gives a solution which has the property that the normal velocity 
w vanishes at the interface; the latter is due to the symmetries of the equations. We 
restrict this solution to R""*"^. We then do the same on the lower halfplane. This 
results in a solution (u2,7r2) for system p.3p that satisfies 

U2 £ i/i(J;Lp(R"+i,R"+i))nLp(J;ff2(K"+i,M"+i))^ 

7r2 e Lp(J;iji(R"+i)), ^2 = on R", 

where, as before, U2 = {v2tW2). We remark that the tangential part of the velocity, 
i.e. V2, may now have a jump at the boundary y = 0. 
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Step 3. To remove the jump in the tangential velocity, we solve the homogeneous 
Stokes problem in the lower halfplane with this jump as Dirichlet datum, that is, 
we solve 

' pidtu — fiiAu + Vtt = 



m 
on 



(3.4) 



pn+l 



div u = 

V = lv2j, W = 

lio(O) = 

where U2 — {v2,W2) is the solution obtained in Step 2. It follows from Proposi- 
tion [3l3] below that system ()3.4p has a unique solution with the regularity properties 
of Theorem 13. II Let (msjTTs) be defined by 

"^''^^ ' \ the solution of dSll) in K'l+K 

Then (u^^tt^) also satisfies the regularity properties stated in Theorem 13. II and we 
have Ivsj = -lv2j and [wg] = on M". 
Step 4. In this step we consider the problem 

pdfU ~ /iAu + Vtt = 



-l^i^yvj ~ 

-2lpdy 



div u 





lf^dy{V2 



f 1^1 = gn 

u(0) = 



gu 



vs)] + [MV:,(ii;2 +W3)1 

- 2lfidy{w2 + Ws)] - ItT2 + TTsj 



m 
in 
on 
on 
on 
in 



pii+i 



pri+l 



(3.5) 

with (v2, 11)2, 7^2) and (v^jW^^tts) the solutions obtained in Steps 2 and 3. Here it 
should be observed that the function on the right hand side of line 3 appearing 
as boundary condition has zero time trace. Problem (|3.5p . which is also of inde- 
pendent interest, will be studied in detail in the next section. It will be shown in 
CoroUarv 14. 21 that it admits a unique solution, denoted here by (U4, n^), which sat- 
isfies the regularity properties stated in Theorem 13.11 

To finish the proof of Theorem l3.1l we set {u, n) ~ (X]^=i where {ui, tt^) 

are the solutions obtained in Step i, with tti := 0. Then (u, tt) satisfies the regularity 
properties stated in the Theorem and it is the unique solution of p.ip . □ 

Remark 3.2. We refer to the recent paper by Bothe and Priiss 9J for results 
related to Theorem 3.1 for the more general and involved situation of a generalized 
Newtonian fluid. 



Let us now consider the problem 

' pdtu — p,Au + Vtt = 
div u — 

U = Uh 

u(0) = 

and prove the result that was used in Step 3 above 



m 
in 
on 
in 



pn+l 



(3.6) 
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Proposition 3.3. Let 1 < p < oo and assume that pj and fj,j are positive constants, 
j = 1,2, and set J = [0, a] . Then problem p.6p admits a unique solution (u, tt) with 

u e oJ?p(J;ip(M"+\R"+')) n Lp{J; H^{W'+\K''+^)), n e Lp(J;iij;(R"+i)) 

if and only if the data Ub = {vt, wt) satisfy the following regularity assumptions 

(a) vb e oWl-^/^P{J;Lp{W\R^))nLp{J:Wp-'^''iW\W')), 

(b) Wb e oH\j;Wp'/''{W'))nLpiJ;W^"'^''{R^)). 



Proof, (i) Assume for a moment that we have a solution in the proper regularity 
class even on the half- line J = K+. Then we may employ the Laplace transform 
in t and the Fourier transform in the tangential variables x £ M", to obtain the 
following boundary value problem for a system of ordinary differential equations 
on M : 

' oj^v - fid^v + i^TT = 0, yj^O, 
bj^w — fJ,dyW + dyTT = 0, y 7^ 0, 
[im + dy^ = 0, y^O, 
■0(0) = Vb, w(0) = Wb. 
Here we have set uJj = pjX + MjlCPi j — li 2, and 



z;,(A,e,y) = (27r)-«/2 



-Ate^*(^l«)w(i, X, y) dx dt, (-l^y > 0. 



This system of equations is easily solved to the result 



V2 




0,2 






11)2 







7r2 






P2A 



for y > 0, and 







ai 






Wl 


















(3.7) 



(3.8) 



for y < 0. Here G K." and have to be determined by the boundary conditions 
v{0) = Vb and w{0) — Wb- We have 



and 



02 — i^a2 = Vb = ai — iS^ai 



UJ2 OJl 

where (a|6) := a^^^ for a, 5 G C". This yields 
Oj = -Ob + j = 1,2, 

+ 



OL2 



ai 



P2A|e| 
PiAiei 



(VA*2(iCI'^6) - W2Wb), 
(V7il(*^l*h) + WiWb). 



(3.9) 
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(ii) By parabolic theory, the velocity u has the correct regularity provided the 
pressure gradient is in Lp, and provided 

Ub e oW^p"'/'^(J;ip(K",K"+i))nLj,(J;Wp2-i/P(R",R"+i)), 

see for instance Denk, Hieber and Priiss [14j . In particular this regularity of Ub is 
necessary. Note that the embedding 

oHliJ; W^p-i/P(R")) n Lp{J; W^-^/PiR'^ j) ^ oW^^'/'^(J; Lp(M")) (3.10) 

is valid. This follows from the fact that Wp"^^^(R") ^ Wp"^/''(R") by a similar 
argument as in the proof of 29, Lemma 6.3] where we set Au := (1 — A)m. 

(iii) We will now introduce some operators that will play a crucial role in our 
analysis. We set G := dt in X := Lp{J; Lp(R")) with domain 

D{G)=oHl{J;Lp{R^)). 

Then it is well-known that G is closed, invertible and sectorial with angle 7r/2, and 
— G is the generator of a Cg-semigroup of contractions in Lp(K"). Moreover, G 
admits an i7°°-calculus in X with i/°°-angle 7r/2 as well; see e.g. [23]. The symbol 
of G is A, the time covariable. 

Next we set := - A, the Laplacian in Lp(M") with domain D(Ai) = i?p(M"). 
It is also well-known that _D„ is closed and sectorial with angle 0, and it admits a 
bounded i/°°-calculus which is even 7?.-bounded with 7?.i/°°-angle 0; see e.g. |15) . 
These results also hold for the canonical extension of Z?„ to X, and also for the 

1/2 1/2 

fractional power Dn of Z3„. Note that the domain of _D„ is 

The symbol of Z3„ is that of d}/^ is given by |^|, where ^ means the covariable 
of X. By the Dore-Venni theorem for sums of commuting sectorial operators, cf. 
[16|,l33j. we see that the parabolic operators Lj := pjG + fijDn with natural domain 

D{L,) - D(G) n D(A.) = oHliJ; ip(M")) n Lp( J; H^iW')) 

are closed, invertible and sectorial with angle 7r/2. Moreover, Lj also admits a 
bounded 77°°-calculus in X with iJ°°-angle 7r/2; cf. e.g. j30]. The same results are 

1 /2 

valid for the operators Fj = Lj , their i7°°-angle is 7r/4, and their domains are 
D(Fj) = D{G'^')nD{Di/^)=oHl/^{J;LpiW'))nLp{J;H;,{R^)). 

The symbol of Lj is PjX + /ijICP and that of Fj is given by ^ pjX + jij 

Let R denote the Riesz operator with symbol C = ^/|^|. It follows from the Mikhlin- 
Hormander theorem that i? is a bounded linear operator on W^(R"), and hence 
also on Lp{J; (R") by canonical extension. 

(iv) Let P2 = P2Xa2- Then the transform of the pressure 7:2 in R"^^ is given 
by e~'^'^/32- The pressure gradient will be in Lp provided the inverse transform of 

132 is in the space Lp{J; Vi^p~^^^(R")). In fact, e^'^'^ is the symbol of the Poisson 

1/2 

semigroup P(-) in Lp(IR"), and the negative generator of P(-) is DJ . Then the 
second part of ({211)) shows that dI/^P{-)(32 e Lp(]R+; Lp(R")) if and only /?2 G 
li'p"^^^(M"). This result extends canonically to Lp( J; Lp(E++^)). 
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Therefore, let us look more closely at /32- We easily obtain 

P2 = P2|^Wh + + M2|$|)(wb - «l«6)), 

where C = HMV We recall that b]!"^ := T-^{\i\T ■) : Wp"(R") W^-^i^"") is an 
isomorphism. 

With the operators introduced above, 62, the inverse transform of /?2, can be 
represented by 

62 = piGD-^I^Wi, + (V7^^^2 + M2^y' )(wfc - ili^li'b)) =: 621 + ^22. 

Due to (jXTOll and oM^p"^/^^( J; n ip( J; W^p"^/^(R")) = £'^,(2 - l/p,p), 

the second term 622 is in 

which embeds into Lp( J; wl^^l^{W^)). Here we used the notation 

Dp^(B,v) = D(F_,-))e,p, (1 + = (D(F,), D(F/))e^p, Q e (0, 1). 

Thus it remains to look at the first term 621 = p2GDn ^^'^Wb- Since 

is bounded and invertible, we see that the condition e qH^{J]Wp ^^^{W^)) is 

necessary and sufficient for 621 G Lp( J; T4^p^^^'°(R")). Of course, similar arguments 
apply for the lower half-plane. □ 

4. The Dirichlet-to-Neumann Operator for the Stokes Equation 

The main ingredient in analyzing problem p.ip with (/, fd, uq) = (0, 0, 0) and 
9v{0) = is the Dirichlet-to-Neumann operator. It is defined as follows. Let {u,tt) 
be the solution of the Stokes problem (|3.6p with Dirichlct boundary condition 
on R", see Proposition 13.31 We then define the Dirichlet-to-Neumann operator by 
means of 

{VN)ui, = -|5(u,^)]e„+i = -|/i(Vu+ (Vu)T)]e„+i -I- He^+i. (4.1) 

For this purpose it is convenient to split u into u = (w, w) as before, and Ub into 
Ub = {vb,Wb)- Then we obtain 

{VN)ub = {-ItJidyv} - [^LVM^~2[^l^yw\ + H). (4.2) 

We will now formulate and prove the main result of this section. 

Theorem 4.1. The Dirichlet-to-Neumann operator VN for the Stokes problem is 
an isomorphism from the Dirichlet space Ub = (vb,Wb) with 

wb e offp(J; vi^i/f (R")) n Lp{j- %2-i/p(r")) 

onto the Neumann space g — {gv,gw) with 

9v e oW^p/'"'/'P(J;ip(M",K"))nLp(J;W^pi-i/P(R",R")), 
g^eLp(J;%i-i/f(R")). 
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Proof, (i) Let (wi , i&i , tti ) and {v2,W2,tt2) be as in ()3.7|) - p.8|) . We may now compute 
the symbol of the Dh'ichlct-to-Neumann operator to the result 

uji^/jr[ai + uj2y/Jl2a2 - ("iMi + c^2^J-2)\S,\^S, - IfJ-ji^m 
2i{p2a2 - AiiailO + 2(a2M2 - aiMi)ICP + A(a2P2 - "iPi) 



iVN)ut 



where the functions aj and aj are given in (|3.9[) . Simple algebraic manipulations 
then yield the following symbol 



where C — and 

a = ^/JI^uJi + ^/JI^uj2, /? = (y"i + M2)|^|, 



7 = {VJj^^2 - VMi^i) 



(4.3) 



(4.4) 



Next we want to compute the inverse of the Dirichlet-to-Neumann operator. Thus 
we have to solve the equation (T>N)uii = g. As before we use the decomposition 
Ub = (wfcjWh) and g = {gv^gw)- Then in transformed variables we have the system 

avb + /3C(CI*b) + ilCwb = 9v, 
-il{C,\vb) + (a + 5)wb = gw 

This yields 

Vb = a-\g, - C(/3(C|«6) + ^7Wfc)]■ (4-5) 

(ii) This last equation shows that it is sufficient to determine (wblC) and Wb- If the 
inverses of /?(t)f,|C) and ^lib belong to the class of g^^ then Vb is uniquely determined 
and has the claimed regularity. Indeed, a is the symbol of 

F VMli^i + Va^F2, D(F) - oHy\j;L^{W^))r\L.p{J-Hl{W)), 

which is a bounded invertible operator from its domain into Lp{ J; Lp(W^)), and 
hence also from Dp{2 — l/p,p) into Df{1 — l/p,p). Here we note that 

Df{0,p) = Dp^ie^p) = oWl/\j;Lp{R")) n Lp{J;W^(R")), 

for 9 € (0,2), 6* 7^ 1. Therefore, F'^g^ belongs to Df{2 - l/p,p) if and only if 

1 /2 

gv € Dp{l—l/p,p). Next we note that 7 is the symbol of ^/]I2F2 — ^fJi\F\ — \[t\F)n 
which is bounded from £'f(2 — '^/p,p) to -0^(1 — l/p,p), and P is the symbol of 

1 /2 

(/ii + fi2)Dn which has the same mapping properties. 

(iii) It remains to show that Wb and {R\vb) belong to Dp{2 — l/p,p). For lib and 
{C\vb) we have the equations 

{a + f3){C\vb) + ilWb = {C\9v), 

-il{C\ib) + {a + S)wb = gw 

since |C| = 1- Solving this 2-D system we obtain 

Wb = ■nir'^[ij{C\gy) + {a + P)gw], 



{(\vb)^m, ^[{a + S){(\gy) -ijgu 



(4.6) 



where 



m 



{a + P)ia + S) - 
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(4.7) 



Since 5 = (3 + {pi + p2)\/\£,\ we obtain the following relation for m 

m={a + I3)[{pi + p2)4 + 4(— + — )^'] (" + 
where 771 = ./piuJi + p2\i\ and 772 = ^/Ji2'^2 + This yields 

(a + p)n n 

(a + pjn n a + p 

We define the operators Tj by means of their symbols rjj, i.e. 

Ti := VmFi + /i2i?y', T2 := + piDlJ\ 0(7,) = D(F,) = D(F) 

Then by the Dore-Venni theorem operators, the operators Tj with domains D(Tj) = 
D{Fj) — D{F) are invertible, sectorial with angle 7r/4. Moreover, they admit an 
_ff°°-calculus with iJ°°-angle 7r/4, see for instance [3D]. The harmonic mean T of 
Ti and T2, i.e. 

T := 2TiT2(Ti + Tz)"! = 2(Tf ^ + T:^^)-^ 

enjoys the same properties, as another application of the Dore-Venni theorem shows. 
The symbol of T is given by 77 := 2771772/(771 +772). 

— 1/2 

Next we consider the operator GDn with domain 

D(Gi?-i/2) = {h^ niD]/^) : D-^/'h e V{G)) 

= oi?p( J; H;\W')) n Lp( J; Lp(K")) 

The inclusion from left to right in the last equality is obvious. The converse can be 
seen as follows. Let h e f)Hl{J; H-^{W'))C^Lp{J; Lp{W^)) and define g := Dn^^^h. 
Then 

g e oi/p(J;ip(M")) nLp(J;iji(M")) ^ Lp{J;HliR")), 

and Dl/^g = ol/^g = h <E Lp[J] LpiW")), which implies that h e TZ{d]J'^) and 

g = Dn^^^h = Dn^^^h G V{G). The operator GDn^^"^ is closed, sectorial and 
admits a bounded _ff°°-calculus with i?°°-angle 7r/2 on X = Lp{J; Lp(]R")); see for 
instance [22l Corollary 2.2]. Its symbol is given by A/|^|. 
Finally, we consider the operator 

N:=ipi+p2)GD-'^^ + 2T, (4.8) 

with domain 

D{N) = D(Gi^-i/2) n D(y) ^ n Lp(J; i7pi(E")); 

recall (|3.10p . By the Dore-Venni theorem N is closed, invertible, and by [30] admits 
a bounded i/°°-calculus as well, with i/°°-angle 7r/2. Its symbol is 77. 

The operator with symbol 7 is then given by T2 — T1, and the operator with symbol 
a + f3 by Ti + T2. For the inverse transforms and {R\vi,) of and (Cl^'b) we 
then obtain the representations 

Wb - N-^[{T2 - Ti){Ti + T2)-^i{R\g^) + <?^] 
(i?|z;fa) = (Ti -I- r2)-'(pi + P2)Gi?-i/2iV-i(fi|5.) (4.9) 
+ N-\R\g,) - {T2 - Ti){Ti + T2)-^N-Hg^ . 
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We note that N ^ has the following mapping properties 

iV-i :ip(J;ip(M"))^ oi?p(J;i?p"'(»"))nLp(J;i7i(M")) -^Lp(J;Lp(R")), 
: Lp(J;iji(E")) ^ oHI{J; Lp{W')) D Lp{J; H^{W')) ^ Lp( J; Lp(R")). 
Therefore by three-fold real interpolation 



Moreover, N-^ maps oM^y^^^^^^l^; ip(IR")) into 



(4.11) 



Next we note that the operators Tj{Ti +T2) ^ are bounded in Di?(l — l/p,p), as is 
the Riesz transform R, and the assertion for Wf, follows now from (|4.9p - (|4.10p and 



The assertions for {R\vb) follow readily from p.lOp and (|4.9p - (|4.1ip 



□ 



We can now formulate our second main result of this section concerning the 
solvability of the problem 



(4.12) 



pdtu — /iAw + Vtt 


= 


in 




div u 


= 


in 


M"+i 




= 9v 


on 


M" 




= 9w 


on 


M" 


I"] 


= 


on 


M" 


u{0) 


= 


in 





Corollary 4.2. Let 1 < p < cx) and assume that pj and pj are positive constants, 
j — 1,2, and set J — [0,a\. Then (|4.12p admits a unique solution (u,7r) with 

u e oifp(J;ip(M"+\M"+i)) nLp(J;i/2(Mn+i^M"+i)), TT e Lp(J;ijpi(M"+i)) 

i/ and only if g = igv,9w) satisfies the following regularity assumptions 

(a) 5. e oW'y'-'/'^(J;^p(K",M"))nLp(J;W^p^-^/''(E",M")), 

(b) g^eLp(J;M^p"'/''(M")). 



Proof. Let u;, := {vb,Wb) :— {T>N)^^{gy,gw), and let (MjTt) be the solution of (|3.6p . 
Thanks to Theorem 14.11 and Proposition 13. 3[ {u,Tr) satisfies the regularity assertion 
of the Corollary, and it is the unique solution of (|4.12D due to the definition of 
VN. □ 

Remark 4.3. The representation formulas in p.Tp ^ p.Sp have also been derived 
and used by other authors, see for instance |11[I34| . However, the optimal regularity 
results in Theorem 13. li Proposition 13. 3[ Theorem 14. 1[ and Corollary 14.21 are new. 
Moreover, the computations and arguments leading to these results are shorter 
than in [TT] (which only deals with the case p — 2) and in [33|. We should mention, 
however, that these authors consider more general domains. 
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5. The Linearized Two-Phase Stokes Problem with Free boundary 
In this section we consider the full linearized problem 



pdtu — fiAu + Vtt 


-/ 


in 




div u 


~ fd 


in 


M"+i 


-Ifidyvj - IfiV^wj 


= 9v 


on 


M" 


< -2li^dywj + Itt] - a Ah 


= 9w 


on 


R" 


M 


= 


on 


R" 


dth — jw 


= gh 


on 


R" 


u{0) = uo, h{0) 


= ho- 







We are interested in the same regularity classes for u and tt as before. Then the 
equation for the height function h lives in the trace space M^p~^^^^( J; Lp(R")) n 
Lp{J; Wp~^^^{W^)), hence the natural space for h is given by 

h e Wp2-i/2p( J; Lp(R")) n Hl{J; VFp2-i/P(R")) n Lp( J; W"3-i/p(r")). 

Our next theorem states that problem (j5.ip admits maximal regularity, in particular 
defines an isomorphism between the solution space and the space of data. 

Theorem 5.1. Let 1 < p < oo be fixed, p ^ 3/2,3, and assume that pj and fij are 
positive constants for j — 1,2, and set J — [0, a] . Then the Stokes problem with free 
boundary (|5.ip admits a unique solution (u, vr, h) with regularity 

u e Hl{J- Lp(R"+i,R"+i)) n Lp(J;i/2(R«+i,R"+i)), 

7reLp(J;ij^i(R"+i)), 

M G Wl/^-^'^P{J- Lp(R")) n Lp{J- W^;-i/P(M»)), 

h e W"2-i/2p( J. 2,p(R")) n Hl(j- vk2-i/p(k«)) n Lp{j- wl-^'p{W')) 

if and only if the data (/, fd, g, gh, uo, ho) satisfy the following regularity and com- 
patibility conditions: 

(a) /eLp(J;Lp(R"+i,R"+i)), 

(b) fd e H^{J;H-HR''+')) n Lp(J;i?i(R"+i)), 

(c) g = {gv,g^o) e l</'"'/'P(J;Lp(R",R"+i)) nLp(J;<-'/^(R",R"+i)), 

(d) gh e Wp-'^^^iJ; ip(R")) n Lp( J; Wp-'^" {W')), 

(e) uo e Typ"'/''(R"+i,R"+i), ho £ Wp'^^'^iW), 

(f) divuo = fd{0) in R"+i and |uo] ^ on R" if p > 3/2, 

(g) -Ipdyvoj - IfiV.wo] = 5.(0) on R" if p > 3. 

The solution map [{f , fd, g, gh,uo,ho) i— > {u,ii,h)\ is continuous between the corre- 
sponding spaces. 

Proof. Similarly as in the proof of Thereom 13.11 we will reduce system (jS.ip to the 
case where (/, fd,g, uq, ho) = (0, 0, 0, 0, 0) and 5^(0) = 0. The Neumann-to-Dirichlet 
operator will once again play an essential role in order to treat the resulting reduced 
problem. 
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(i) Let 



-ire-(i+A.)*_ 



-2(1+D„)tl 



where uq = {vo,Wq) and 7 : 
the following properties 



is the trace operator. The function hi has 



/ii(0) = /io, at/ii(0) =5h(0) + 7wo, 



(5.3) 



see [ini Lemma 6.4] for a proof of a similar result. Let then {ui, tti) be the solution 
of problem (13. ip . with g^, replaced by 17^1 + crA/ii. It follows from Theorem (23 the 
assumptions on g = {gv,gw), and from the first line in (j5.3p that (ui,7ri) satisfies 
the regularity properties stated in Theorem 15.11 

(ii) Next we consider the reduced problem 



pdfU — fiAu + Vtt 


= 


in 




divM 


= 


in 




-Ifidyvj - IfiV^wj 


= 


on 


M" 


2lfidywj + Itt] - a Ah 


= 


on 


M" 




= 


on 


M" 


dth — 


= .9/1 


on 




u(0) = 0, h{0) 


= 0, 







(5.4) 



with g/i := .g/i— (c)t/ii — 7W1), where ui = (wijiyi) is the solution obtained instep (i). 
We conclude from (|5.3p and the regularity properties of 7101 that 

h e oW^-'^'PiJ; Lp(M")) n Lp(J; ^^-^/f (M")). (5.5) 

Suppose that problem (|5.4p admits a solution (m2, 7r2, /12) with the regularity prop- 
erties stated in (j5.2p . One readily verifies that (m, tt, /i) := (ui + U2, tti + 7r2, hi + /12) 
is a solution of problem (|5.ip in the regularity class of (|5.2p . 

(iii) It thus remains to show that the reduced problem (|5.4p admits a unique so- 
lution {u, TT, h) in the regularity class stated in Theorem 15.11 We note that once 
h has been determined, Corollary 14.21 yields the corresponding pair (w,7r) in prob- 
lem 

To determine h we extract the boundary symbol for this problem as follows. Ap- 
plying the Neumann-to-Dirichlet operator {VN)^^ to {gv,9w) — {0,crDnh) yields 
ju — Uh, the trace of u. According to (|4.7p . the tranform of the normal component 
jw = Wb of Ufc is giyen by 



Wb 



-h. 



{pi +P2)A/|C| +4?7i 772/(771 -^772) 
Let us now consider the equation dth — = gh- Inserting this expression for Wb 
into the transformed equation Xh — Wb — gn results in s(A, |f |)ft. — gn where the 
boundary symbol s(A, |^|) is giyen by 



s(A,|ei) = A + 



(Pl+P2)A/|e| +4771 772/(771 +??2)' 



(5.6) 
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The operator corresponding to this symbol is 

S = G + (jDnN-^, (5.7) 

where the meaning of the operators G, £>„ and N is as in Section 4. S has the 
following mapping properties: 

S : oHl+\j- K;{R-)) n oH;{J; K;+\R-)) -> oH'-iJ; K^iR^)), (5.8) 

where K S {H, W}. In order to find h we need to solve the equation Sh = gh, that 
is, we need to show that S is invcrtible in appropriate function spaces. 

All operators in the definition of S commute, and admit an i/°°-calculus. The 
H°°-a.ngle of Dn is zero, that of iV is 7r/2 and that of G is 7r/2 as well. Thus we 
can a-priori not guarantee that the sum of the power-angles of the single operators 
in 5* is strictly less than tt, and the Dore-Venni approach is therefore not directly 
applicable. We will instead apply a result of Kalton and Weis [Ml Theorem 4.4]. 

For this purpose note that for complex numbers Wj with argWj e [0,7r/2), we 
have a,rg {W1W2) / {wi + W2) ~ arg(l/z«i -f- l/w2)^^ G [0, 7r/2) as well. This implies 
that s(A, 1^1) has strictly positive real part for each A in the closed right halfplane 
and for each e K", (A, £_) ^ (0, 0), hence s(A, |^|) does not vanish for such A and ^. 

We write s(A, |f |) in the following way: 

s(A,t) = A + CTTfc(z), z = A/r^, A G C, T G C\ {0}, (5.9) 

where 

k{z) = [(pi+p2)2 + 4( 



The asymptotics of k{z) are given by 

fc(0) = — — -, zk(z) for z eC\R- with Izl 00. 

2(/ii +P2) P1+P2 

This shows that for any i3 G [0, tt) there is a constant C — C(i?) > such that 

|fc(z)|<^^, zGS^. 

Hence we see that 

|s(A, 1^1)1 <C(|A| + 1^1), ReA>0, ^gR", 

is valid for some constant C > 0. Next we are going to prove that for each Aq > 
there are 77 > 0, c > such that 

|s(A,t)| >c[|A|-f |r|], for all A G E,/2+^, |A| > Ao, T G S^. (5.10) 

This can be seen as follows: since Refc(z) > for Rez > 0, by continuity of the 
modulus and argument we obtain an estimate of the form 

\s{X, t)\ > co[|A| + |rp(z)|] > c[|A| + |t|], A G E,/2+^, r G 

provided \z\ < M , with some 77 > and c > depending on M, but not on A and r. 
On the other hand, for m > fixed we consider the case with |A| > m\T\, \z\ > M. 
We then have 

\s{X, r)| > |A| - a|r||fc(z)| > i[|A| + m|r|] - aC\r\/{l + M) > c[\\\ + |r|], 

provided m > 2(7(7/(1 + M), and then by extension 

|s(A, t)| > c[|A| + |r|], A G S,/2-f,„ r G |A| > m|r|, \z\ > M, 
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provided > and c > are sufficiently small. One easily sees that the intersection 
point of the curves y — Mx^ and y = mx in M? has distance d — {m/M)y/l + m? 
from the origin. By choosing M large enough so that d < Ao, (IS.lOp follows by 
combining the two estimates. 

By means of the 7?,-boundedness of the functional calculus for Z)„ in /^^(M"), cf. 
Desch, Hieber and Priiss [TS], we see that 

{X + D]/^)s-\\^dU^) 

is of class H°° and 7?.-bounded on S7r/2+jj \ -Bao(O)- The operator-valued H°°- 
calculus for G = dt on oH'p{J; Kp{R^ j), cf. Hieber and Priiss [53], implies bound- 
edness of 

{G + Di/^)s-\G,Dl,/^) in oH;{J; K^iR^)). 
This shows that s^^{G, dI/^) has the following mapping properties: 
s-\G,D'J') : oH;{J;K;{W')) ^ oH'p^\J; K;{Wl)noH;{J; K;+\R-)). (5.11) 

We conclude that S is invertible and that S^^ = s^^{G, dI/'^). Choosing r = and 
s = 2 - 1/p and K ^ W in ((5TT|) yields 

: Lp( J; %2-i/p(M")) ^ oi?p( J; n Lp{J; M^p^-i/f (M")). (5.12) 

Moreover, we also obtain from (|5.1ip 

: Lp(J;Lp(R")) offp( J; ip(K")) 
: oHl{J;Lp{R'')) ^ oi^p( J; ip(M")). 
Interpolating with the real method (• , ■)i-i/p^p then yields 

: oWl-^^P{J; Lp(R")) oW^p"'/^( J; ip(M")). (5.13) 

(|5.12p - (|5.13p shows that the equation Sh — gh has for each gh satisfying (|5.5p a 
unique solution /i in the regularity class (|5.2p . 

(iv) Since the function is now known we can use CoroUarv 14.21 to determine the 
pair (u, tt) in problem (|5.4p . For this we note that 

see [291 Lemma 6.2] for a proof. This shows that the function h determined in step 
(iii) satisfies 

Ah e oW^/^-^/^P{J; Lp(K")) n Lp{J; ^^^^-^/^(M")) 

and Corollary 14. 21 vields a solution {u,tt) in the regularity class (|5.2p . 

(v) Steps (i)-(iv) render a solution (w, tt, /i) for problem 15. ll that satisfies the regu- 
larity properties asserted in the Theorem. It follows from step (iv) and from The- 
orem [ST] that problem (|5.4p with {f,fd,g,gh,uo,ho) = (0,0,0,0,0,0) has only the 
trivial solution, and this gives uniqueness. The proof of Theorem 5.1 is now com- 
plete. □ 

Remark 5.2. Further mapping properties of the symbol s(A, r) and the associated 
operator S have been derived in [31] . In particular, we have investigated the sin- 
gularities and zeros of the boundary symbol s, and we have studied the mapping 
properties of S in case of low and high frequencies, respectively. 
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6. The nonlinear problem 

In this section we derive estimates for the nonhnear mappings occurring on 
the right hand side of (|2.9p . In order to facihtate this task, we first introduce 
some notation, and then study the mapping properties of the nonhnear functions 
appearing on the right hand sight of equation (j2.9|) . In the following we set 



(6.1) 



Ei(a) :={uei7i(J;ip(M"+i,M"+i))nLp(J;il2(R"+i,M«+i)) : H = 0}, 
E2(a) :=Lp(J;ijpi(M"+i)), 

E3(a) I^;/2-i/2P( J; n Lp(J; (M")), 

E4(a) := W^p2-i/2f(J;Lp(R"))ni7pi(J;I^2-i/p(j^n)) 

n T4^pi/2-i/2p( J. H^(^-n-^) n Lp{J; iy3-i/f (R")), 
E(a) := {{u,TT,q,h) E Ei(a) x E2(a) x £3(0) x E4(a) : [tt] = g}. 
The space E(a) is given the natural norm 

\\{u,Tr,q, h)\\E{a) = ||"||Ei(a) + lk||E2(a) + ll9l|E3(a) + WHEiia) 

which turns it into a Banach space. We remind that ¥.2(0.) is equipped with the 
norm = (EJ^^ l|5,-r,^(,,,^(«„,.)//^ for n : ^ M. 

In addition, we define 

Fi(a) :=Lp(J;Lp(R"+i,M"+i)), 

F2(a) iIpi(J;ii-i(R"+i)) n Lp( J; ili(R"+i)), 

F3(a) W^pi/2-i/2p^j.^^(j^„^jg«+i))^^^(j.p^^i-i/p^jg«^jjn+i))^ (g 2) 

F4(a) W^'^/^P{J; ip(R")) n Lp( J; W^2-i/P(Rn))^ 
F(a) := Fi(a) x F2(a) x W-sia) x F4(a). 

The generic elements of F(a) are the functions {f,fd,g,gh)- 

We list some properties of the function spaces introduced above that will be 
used in the sequel. In the following we say that a function space is a multiplication 
algebra if it is a Banach algebra under the operation of multiplication. 

Lemma 6.1. Suppose p > n + 3 and let J = [0, a]. Then 

(a) E3(a) andW^i^a) are multiplication algebras. 

(b) Ei(a) ^ C(J;B?7Ci(K"+\K"+^))nC(J;Si7C(R"+i,R"+i)). 

(c) E3(a) ^ C(J;BUC(R")). 

(d) E4(a) SCi(J;BCi(R"))nBC'(J;SC2(M")). 

(e) W^~'^^P{J; Lp(E")) n H^iJ; W^-'^^iW^)) n Lp{J; W^'^^iW^)) ^ Ei{a). 

Proof, (a) The assertion that E3(a) and F4(a) are multiplication algebras can be 
shown as in the proof of |:29., Lemma 6.6(ii)]. 

(b) It follows from m Theorem III.4.10.2] that Ei(a) C( J; Wp^"^/''(R"+\ R"+i)) 
and this implies the first inclusion, thanks to Sobolev's embedding theorem. The 
second assertion follows from the fact that u is continuous across y — 0. 

(c) This follows from [191 Remark 5.3(d)] and Sobolev's embedding theorem. 
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(d) We infer from [21 Theorem III.4.10.2] that 

and the inclusion £4(0) C(J; i3C^(M")) then foUows from Sobolev's embedding 
theorem. In addition, we conclude from [^Hl Remark 5.3(d)] and Sobolev's embed- 
ding theorem that 

Wp^-2/P(J;Lp(K"))nLp(J;Ty2-i/P(Rn)) ^ BC{J; BCHM.'')) 
and this implies that E4(a) ^ BC^{J; BC^W"-)). 

(e) This follows from (jSH)) . □ 
Let 

iV(u, TT, g, /i) tt, h),Fd(u, h), G{u, q, h), H{u, h)) (6.3) 

for (u,7r, g, /i) g E(a), where as before u — {v,w), F = [F^,Fw) and G — {Gv,Gw). 
We show that the mapping N is real analytic. 

Proposition 6.2. Suppose p > n + 3. Then 

N e G'^{E{a) ,¥{a)) and A^(0) = 0, L»Ar(0) = 0, (6.4) 

where DN denotes the Frechet derivative of N . In addition we have 

DN{u, TT, q, h) e £(oIE(a), oF(a)) for any {u, tt, q, h) G E(a). 

Proof. We first note that the mapping [iu, tt, q, h) 1— s- A^(u, tt, g, /i)] is polynomial. 
It thus suffices to verify that N : E(a) — > E(a) is well-defined and continuous. 

(i) We first consider the term F{u, tt, h), and observe that it contains the expressions 
Wh,Ah and dth. Without changing notation we here consider the extension of h 
from R" to R"+i defined by h{t, x, y) = h{t, x) for t e J and (a;, y) e M" x M. With 
this interpretation we clearly have 

l|5/i||oo„/xR"+i - l|5/i||co„7xR", /iGE(a), 5e {5j,A,5t}, (6.5) 

where || • \\oo,u denotes the sup-norm for the set [/ C J x M"^-'^. Next we note that 

BC( J; BC(R"+i)) • Lp(J; ij,(M"+i)) ip( J; Lp(M"+i)), 

BC{J\ BC(R"+^)) • BC( J; BC(R"+^)) BC(J; BC(R"+^)), 

that is, multiplication is continuous and bilinear in the indicated function spaces. 
We can now conclude from (I6.5I)-(I6.6I) and Lemma [O] that 



F e C"(Ei(a) X E2(a) x E4(a),Fi(a)), F(0) = 0, DF{Q) = 0. 

(ii) We will now consider the nonlinear function Fd{u, h) = {yh\dyv). Since h does 
not depend on y we have 

Fa{u, h) = {\Ih\dyu) = dy{\Ih\u). (6.7) 

Observing that 

BC\J-BG{W'+^)) ■ Hl{J;L.p{«^+')) ^ H^^iJ; Lp{W^+')), 
BC{J;BC\W'+'}) ■ Lp(J;i/i(R"+i)) ^ Lp( J; ffpi(R"+i)), 

and 

dy e £(iIpi(J;Lp(R"+i)),i/i(J;iJ-i(R"+i))) 
n /:(Lp( J; iJpi(R"+i)), Lp( J, ip(R"+^))) , 
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we infer from Lemma l6. If d) that 

FdGC'^(Ei(a)xE4(a),F2(a)), Fd(0) = 0, i?Fd(0) = 0. 

(iii) We remind that 

l/ia. •] e /:(i?pi(J;Lp(M"+i))nLp(J;i72(M»+i)),E3(a)) (6.8) 

where l/ic^iw] denotes the jump of the quantity fidiU with u a generic function from 
R"+^ R, and where di ~ dx^ for i = 1, . . . , n and dn+i — dy. 
The mapping G(u, q, ft,) is made up of terms of the form 

{^i^^Ukl^jh, IfidiUkjdjhdih, qdjh, Ahdjh, G,^{h), Gt^{h)djh 

where Uk denotes the fc-th component of a function u G Ei(a). From (j6.8p and the 
fact that E3 (a) is a muhiphcation algebra foUows that the mappings 

{u,h) ^ l^idiUk\djh, IfidiUkjdjhdih : Ei(a) x £4(0) E3(a), 

{q, h) ^ qd-jh : £3(0) x £4(0) ^ £3(0), ft h-> Aft^^ft : £4(0) £3(0) 

are multihnear and continuous, and hence real analytic. The fact that £3(0) is an 
algebra additionally implies that the mapping [ft i— > GK(ft)] : £4(0) IE3(a) is 
analytic. In summary we conclude that 

G £ C"(IEi(a) X E3(a) x £4(0), £3(0)), G(0) = 0, DG{0) = 0. 

(iv) We infer from 7 e C{H^{J] Lp{R"+^)) n Lp(J; ij2(K"+i)), Fa(a)) and the fact 
that F4(a) is an algebra that the mapping [(u, ft) ^-^ (Vft|7u)] : Ei(a) x £4(0) — > 
F4(a) is bilinear and continuous. This immediately yields 

H e G"(Ei(a) X E4(a),F4(a)), H{0) ^ 0, DH{0) = 0. 

(v) As the terms of N are made up of products of u, tt, q, ft and derivatives thereof, 
one easily verifies that 

DN{u,Tr,q,h)[u,Tr,q,h] G a¥{a) whenever (it, 7r,(j, ft) e E(a), (m, 7f,q, ft) £ oJE(i)- 

Combining the results obtained in steps (i)-(v) yields the assertions of the propo- 
sition. □ 

We are now ready to prove our main result of this section, yielding existence and 
uniqueness of solutions for the nonlinear problem (j2.9p . 

Theorem 6.3. (Existence of solutions for the nonlinear problem (|2.9p ). 

(a) For every to > there exists a number e = e(to) > such that for all 
initial values 

[uq, fto) e VFp2-2/P(K"+i^R"+i) X W^p3-2/p(M«), |uo] = 0, 
satisfying the compatibility conditions 
|/iI?(uo,fto)t'o - M(^'o|2'(wo,fto)i^o)i^ol = 0, divMo = -Fd(uo, fto), [uq] = (6.9) 
and the smallness condition 

||uo|lvi/2-2/p(j,„ + i^ + ll'*0|lvi^3-2/P(R„) < £, (6.10) 

where Vlu^h) is defined in (|2.6p . the nonlinear problem p.9p admits a 
unique solution (u, tt, [tt], ft) G E(io). 



26 



JAN PRUSS AND GIERI SIMONETT 



(b) The solution has the additional regularity properties 

(u,7r) e C"^((0,to) X M«+\M"+2), |7r],/i e C"^((0,to) x M"). 

/n particular, M = Ute(o io) ('i^J' ^ ^l^)) '"^'^^ analytic manifold. 

Proof. In order to economize our notation we set z := (m, tt, q, ft,) for (u, tt, q, h) £ 
E(a). With this notation, the nonhnear problem (|2.9p can be restated as 

Lz^N{z), KO),/i(0)) = K,/io), (6.11) 



where L denotes the hnear operator on the left-hand side of (|2.9p . and where N is 
defined in (|6.3p . 

It is convenient to first introduce an auxiliary function z* G E(a) which resolves 
the compatibility conditions (j6.9p and the initial conditions in (j6.1ip . and then to 
solve the resulting reduced problem 

Lz ^ N{z + z*) - Lz* =: Ko{z), z e oE{a), (6.12) 

by means of a fixed point argument. 

(i) Suppose that the initial values (uq, Hq) satisfy the (first) compatibility condition 
in (|6.9p . and set 

I^ol := Mvo\T^iuo, ha)iya)j + a{Aho ~ G«(/io)). 

It is then clear that the following compatibility conditions hold: 

-|/i9yt;o]-lAiV,u;o]=a(wo,Nl,/Jo) on R" 

-2lfidyWoj + [tto] - aAfto = G^uo, ho) on R" ^ ' 

where uq = [vq, wq). Next we introduce special functions (0, g*, 5^ ) e IF(a) which 
resolve the necessary compatibility conditions. First we set 

'' 7^+e-*■°"+l£+(t;o|V/^o) in R';+^ 

7^_e"*-°"+lf_(^;o|V/^o) in 

where £± £ C{Wp ^^''(R!^^"'^), Wp "^^^{R"^^)) is an appropriate extension operator 
and TZ± is the restriction operator. Due to (t;o|Vfto) G W^~^^^(R"+^) we obtain 
c* e iJpi(J;Lp(R"+i)) nLp(J;ij2(R"+i)). 

Consequently, 

f*a := dyc* e F2(a) and /rf*(0) = Fd{vo,ho). (6.15) 

Next we set 

g*{t):=e-''-'G{uo,l7Tolho), gl{t) := e-""-* H{uo, h,). (6.16) 

It then follows from (jOS]) and [19, Lemma 8.2] that (0, 5*, g^;) G F(a). ((6l3l) 
and the second and third condition in (|6.9p show that the necessary compatibility 
conditions of Theorem lS.ll are satisfied and we can conclude that the linear problem 

Lz* ^{OJlg*,gl), {u* (0), h* {0)) ^ {uo, ho), (6.17) 

has a unique solution z* G E(a). With the auxiliary function z* now determined, 
we can focus on the reduced equation (|6.12p , which can be converted into the fixed 
point equation 

z^Lo^Koiz), ^GoE(a), (6.18) 



c*{t):={ ^ ' ' (6.14) 
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where Lq denotes the restriction of L to oIE(a). Due to the choice of {f^,g* ,g'^^) we 
have Ko{z) E o¥{a) for any z e oIE(a), and it follows from Proposition 16.21 that 

Ko e C"(oIE(a),oF(a)). 

Consequently, Lq^Kq : oE(a) ^ oIE(a) is well defined and smooth. 

(ii) In the following, to > is a fixed number. We set 

El := {iuo,ho) e W^^'^^PiW'+\W"+^) x W^'^^p(R'') : fuj = 0}, 

and observe that Ei is a Banach space. Given {ua,ho) G Ei let {fd: 9* ^ 9h) be 
defined as in (|6.15p - (|6.16p . It is not difhcult to see that the mapping 

i^* :i?i^F(<o), F*{u^,h„):^{Q,f2,9\9l). 
is (in fact real analytic), and that F*(0) — and DF*{Q) = 0. Hence given 
6 € (0, 1) there exists e = e{S) > such that 

\\F*{uo.ho)\\rito)<S\\iuo,ho)\\E„ (uq^o) e eBb,. (6.19) 

Let G(to) denote the closed subspace of ¥{to) x Ei consisting of all functions 
(/) fd, 9, 9h, uq, ho) satisfying the compatibility conditions of Theorem l5.1l 

Suppose that (ito,/io) G £^Ei satisfies the compatibility conditions (j6.9[) . Then, 
due to (|6.13p and the definition of F* , the mapping 

G* : Ei-^ G{tn), G*{uo,ho) := {F* {uq, ho),Uo, ho), 

is well-defined and ||G'*(mo, /lo)llG(to) — 2|| (uo, /io)||£;i • It then follows from Theo- 
rem lS.ll that ()6.17p has a unique solution z* = z*{uo,ho) which satisfies 

\\z*\\E{to) < Co\\{uo,ho)\\Ei, {uo,ho) e eMe, , (6.20) 
where the constant Co does not depend on (uo,/io). 

(lit) Theorem 15.11 also implies that Lq : oIE(^o) o^ito) is an isomorphism. Let 
then 

M.= \\L^ ||£(oF(to).oE(to))- (6-21) 

We can assume that the number 6 in step (ii) was already chosen sufficiently small 
such that 

'<^^^'^M(^)^- ^'-''^ 

(iv) We shall show that the fixed point equation (|6.18p has for each initial value 
(uo,/io) satisfying ()6.9|) - (l6.10p a unique fixed point z = z{uQ,ho) € eB„E(to)- It 
follows from Proposition 16.21 and (I6.20p that 

\\DN{z + z*)\\c(E(to).¥{to)), \\DKoiz)\\c(oKit„).o¥it„)) <S (6.23) 
for all (mq, ho) satisfying (|6.9p - (|6.10p and all z e eBpE(to)j provided e is chosen small 
enough. From (|6.19p - (|6.23p follows for z, zj e eM^E{to) 

\\Lo\Ko{zi) - Ko{z2))\\„Eito) < MS\\zi - Z2||oE(to) < (l/2)||zi - ^2|loE(to) 

and 

11^0 '^o(2)||oE(to) < M{\\Niz + Z*)\\Eito} + \\F*{uo, ^o)llF(to)) 

< Md{2 + Co)e < e. 

This shows that the mapping Lq^Kq '■ eBpE^j^-) —> eRgE{to) ^ contraction for any 
initial value {uo,ho) satisfying (|6.9p ~ (|6.10p . 
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(v) By the contraction mapping principle L'^^^Kq has a unique fixed point z £ 
£®oE(to) oIE(io) and it follows from (|6.1ip - (|6.12p that z + z* is the (unique) 
solution of the nonlinear problem (j2.9p in E(to)i proving the assertion in part (a) 
of the Theorem. 

(vi) In order to show that (m, tt, q, h) is analytic in space and time we can use the 
same strategy as in [ini Section 8] . Since the proof is similar we will refrain from 
giving all the details, and will rather point out the underlying ideas. 

Let {u,TT,q,h) G lE(to) be the solution of (|2.9p with initial value (uo,ft.o)- Let a E 
(0, to) be fixed and choose ^ > so that (1 + S)a < to. Moreover, let ip he a. smooth 
cut-ofi' function with ip = 1 on [—R,R] for some R > and suppose that (5 > is 
chosen small enough so that 

l + ^{y)Tt>0, 1 + > 0, te[0,a], Te{-S,S), yeR. 

For given parameters (A, z^, r) e (1 — (5, 1 + (5) x M" x {—S, 5) we set 

iux.,i',T, '^\,u,T){t, x,y) : = {u, 7r)(At, x + tv, y{l + ip{y)Tt)), 

iqx,i^,hx,^)it,x) : = {q,h){\t,x + tiy), (6.24) 

zx,u,T ■ = {ux,i^,T,Trx,u,T,qx,u,h\^i,) 

where {t, x, y) E [0, a] x M" x M. Suppose we know that 

[{X,i^,t) >-> ZA,.,r] e C-(A,E(a)) (6.25) 

with A C (1 - (5, 1 + (5) X R" X (-5, 6) a neighborhood of (A, r) = (1, 0, 0). Pick 
{so,xo,yo) & (0,io) x M"^^ and choose a G {so,to)- Without loss of generahty we 
can assume that yo E [—R, R]- Thanks to the embeddings 

Ei(a) ^ C(/,BC(R"+\M"+^)), E3{a),Ei{a) ^ C(/; BC(M")), 

see Lemma l6. 11 we conclude that 

[(A,i.,T) K^iiA,.,r] e C7'"(A,C7(/;BC7(R"+\M"+i)), 

^ (<?„., /lA,.)] e C"(A,C(/;BC(R") x C(/;BC(R")) 

for / = [0,a]. Thus 

[{\,i^,t) ^ u{Xso,xo + sniy,yo{l + Tso)] E C"(A,M"+i), 

[{X,i^) ^ {q,h){Xso,xo + soi^)] E C"(A,R2), 

and this implies that 

ue C"^((0,to) X M"+\R"+i), q,hEC'^{{0,ta)xW). (6.26) 

This in turn with ^^^-^^^ shows that Vtt E C"^((0,io) x R"+\R"+i) as weU, 
and we can now conclude that 

7reC""((0,io) xR"+i), (6.27) 

where the pressure tt is normalized by 7r(t, 0, 0—) = 0, i.e. 

I /o^[(^2:^(^>«2:, sy)|x) +5y7r(t, sx, sy)y]rfs, y < 0- 



(vii) We will now explain the steps needed to establish the crucial property (|6.25p . 
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First we note that there exists a neighborhood A C (1 — (5, 1 + (5) x M" x {—6, S) of 
(1,0,0) such that 



[(A,j.,r) ^ (0,/a,.,.,5l.,<A,.)] e C'^(A,F(a)) 



(6.28) 



where the functions {fd:9*7 9h) are defined in (|6.15p - (|6.16p . In fact, the assertion 
follows immediately from [HI Lemma 8.2] for the functions {g*,g'^). Let us then 



consider the function c* defined in (|6.14p . Let w{t) := e 



'^wq for some function 



wq e Wj 



2-2/p/-|n,„+i 



) and define w\^^^r{t, x, y) for (i, x,y) € / x R"+^ as above, with 



I — [0, a]. Then one verifies as in the proof of [TH Lemma 8.2] that 
wx.^^^r e ffpi(/;Lp(K"+i)) n ip(/;i/2(M"+i)) Xi(/) 

for (A, i^, r) € (1 — 5, 1 + (5) x M" x (—5, (5), and that iua.i/.t solves the parameter- 
dependent parabolic equation 

dtU — Ax^u^rU — 0, u(0) = Wo, 

in R"+^, where ^Ia.i/.t is a parameter-dependent differential operator given by 



A 



a(y) 



Aa"(j/)i 



9y + {.v\Vx) 



{l + a'{y)TtY y \l + a'{y)Tt (1 + a'(?/)rt)3 

for t G [0, a] and y G M, where «(?/) := yip{y). Here we observe that 

A,o,o = A, [(A,z.,t) ^ Ax,,,r] e C"(A,/:(Xi(/),Xo(/)), 

with Xo(/) := Lp(/,Lp(M"+i). As in the proof of [19, Lemma 8.2] it follows from 
the implicit function theorem that there exists a neighborhood Ac(l — (5, l + (5)x 
R" X (-6,6) of (1,0,0) such that 

[(A,i^,r) ^ wx,,,r] e C"(A,Xi(/)). (6.29) 

Applying (|6.29|) separately to wq — £±{vo\/ho), an then applying TZ± yields 

[(A, z., r) ^ c^,,,] e C"(A, i7pi(/; Lp(R"+i)) n Lp(/; i/2(M"+i)). 

It then follows from the definition of that [{X,v,t) fd x u t] ^ (^"(A, F2(a)). 
In a next step one verifies that the function ^ ^ solves the linear parameter- 
dependent problem 



pdtu - Axm,tU + Sa.tTt 


= 


in 




CrU 




in 




1 + rr 




on 




+ rr 


* 


on 


M" 




= 


on 








on 


R" 


m(0) = Mo, /i(0) 


= /lO 







(6.30) 
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where 



(6.33) 



Sa.tTt := A(V:i7r, — — ] dyn), CrU div^v + — — ] dyW, 
l + a'(y)Tt l + a'(y)Tt 

V^h := -{v\Vh). 
We note that 

^1^0,0 = AiA, Buo = V, Ci = div, Vo = 0. 

It is easy to see that the differential operators Ax,^^t, 13\^r, Cr and I?^ depend 
analytically on the parameters (A, r) in the appropriate function spaces. Using 
Thereom lS.ll and the implicit function theorem one shows similarly as in [19, Lemma 
8.3] that there is a neighborhood A C (1 - (5, 1 + (5) x M" x (-6, 6) of (1, 0, 0) such 
that 

[{\,i,,T)^zl,JeC'^iAMa))- (6.31) 
Let z be the solution of (|6.12p obtained in step (v) above. Then one verifies that 
z\,v,T G 2£BpE(jp) for (A, v, t) £ A, with A a sufficiently small neighborhood of 
(1,0,0). Moreover, za.i/.t solves the nonlinear parameter-dependent problem 

Lx^„^rZ ^ K\^^,r{z), zeoE(a), (6.32) 
for (A, z^, r) S A, where L\^^^rZ is defined by the left-hand side of (I6.30p and where 
/ Ai^r(w + wl^,^,7r + 7r*_^^,/i + /i* ,,) \ 
Fd,r{u + w*,^,^ , h + hl^^) - fl^^^^ 
Gr{u + „^ , q + ql^ , h + hl^^) - gl^^ 

The functions i<V, ^d,T and Gr are obtained from F, Fd and G, respectively, by 
replacing terms containing partial derivatives dy and dy in the following way: 

T 1 02 1 o2 ci"{y)Tt 

<Jv^ ^ -. 7T-y OyUJ, 0,,U} ^ — — r7TO,,UJ — -— T:rOyU} 

^ l + a'{y)Tt ^ y {l + a'{ii)Ttf ^ (1 -f- a'(y)rt)3 ^ 

for Lo G {v, w, tt}. Equation (|6.32p can be reformulated as 

^{z,{X,i^,t)) -.^ z-{Lx,,.rr^Kx,,.r{z)--0, zeoHa)- (6.34) 

Here we observe that ^'(z, (1,0,0)) — for the solution z of the fixed point equa- 
tion ((08)) . It follows from ((6T28)) . fOTjl and Proposition [Ql that 

[(z, (A, i^, r)) vl/(z, (A, i^, r))] e C"(oE(a) x A, oE(a)). 

Moreover, it follows from ((OT|) - (p?^ that 

Di^iiz, (1, 0, 0)) = / - D{L-^Ko){z) e Isom(oE(a), oE(a)). 

By the implicit function theorem there exists a neighborhood Ac (1 — 6,1 + 5) x 
W X (-(5, S) of (A, zy, t) = (1, 0, 0) such that 

[(A, 1^, t) ^ ZA,.,r] e C"(A, oE(a)). (6.35) 

Combining (|6.3ip and (|6.35p yields (|6.25p . This completes the proof of Theorem l6.3l 

□ 
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Proof of Theorem 1.1: We first observe that the compatibihty conditions of 
Theorem 1.1 are satisfied if and only if ()6.9|) is satisfied. Next we note that the 
mapping Qho given by (cc, y) := {x,y + ho{x)) defines for each ho g Wp'^^^{M.") 
a C^-diffeomorphism from M."^^ onto f2i(0) with det[£'8ft,o(x, y)] = 1. Its inverse 
is given by &J^^{x, y) :— {x,y — ho{x)). It then follows from the chain rule and the 
transformation rule for integrals that 

where C{h()) := M[l + || V/io||bci(R")], with M an appropriate constant. Con- 
sequently, there exists Eq > such that ||wo|j^2-2/p^j-j^j + ||/io||^3-2/P(-jj„^ < £o 
implies the smallness-condition (|6.10p . Theorem 16.31 then yields a unique solution 
(w, w,7r, [TT],h) £ E(to) which satisfies the additional regularity properties listed in 
part (b) of the theorem. Setting 

{u,q)(t,x,y) = {v,w,Tr){t,x,y - h(t,x)), {t,x,y) € O, 

we then conclude that (m, q) e C"^(0, M"+2) and [q] G C^iM). The regularity prop- 
erties listed in l|1.3P " (|1.4p are implied by Lemma lOT bl-fc). Finally, since n{t, x, y) is 
defined for every {t,x,y) € O, we can conclude that 7r(i, •) g Hp{Q{t)) C UC{n{t)) 
for every i e (0, to)- D 
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